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Abstract. We present hard X-ray and neutron diffraction measurements on the pofar ffuorocarbons 
HCF3 and H3CF under supercriticaf conditions and for a range of molecular densities spanning about 
a factor of ten. The Levesque-Weiss-Reatto inversion scheme has been used to deduce the site-site 
potentials underlying the measured partial pair distribution functions. The orientational correlations 
between adjacent fluorocarbon molecules - which are characterized by quite large dipole moments but 
no tendency to form hydrogen bonds — are small compared to a highly polar system like fluid hydrogen 
chloride. In fact, the orientational correlations in HCF3 and H3CF are found to be nearly as small as 
those of fluid CF4, a fluorocarbon with no dipole moment. 



PACS: 61.20.Qg, 61.20.Ja, 61.25.Em 

1. Introduction 

The understanding of dielectric properties resulting from orientational correlations [1, 2], as well as the 
determination of orientational correlations from diffraction experiments [3, 4], are long standing problems 
in the physics of molecular fluids. The simple fluorocarbons HCF3 and H3CF are very interesting model 
substances in this context, as they posses rather large dipole moments (1.65 • 3.336 • 10~ 30 Cm), in the 
case of trifluoromethane - the same as that of the water molecule) but no tendency to form hydrogen 
bonds [5]. The investigation of the structure of the simple fluorocarbons H3CF and HCF3 thus enables 
the study of the structural effect of the molecular dipole alone. 

Although the properties of these fluorocarbons are interesting, only very limited structural 
information is available so far. HCF3 is discussed as replacement for chlorinated hydrocarbons as 
refrigerant since it it has no ozone damaging effect, it has a shorter atmospheric lifetime and, hence, a 
lower global warming potential and it presents no toxological risk [6, 7]. HCF 3 is discussed for extraction 
applications [8] and it has been shown, that the enatioselectivity of asymmetric catalysis can be controlled 
by the density of the fiuoroform solvent [9]. The crystal structure of HCF3 has been determined by a 
neutron powder diffraction experiment [10] and the molecular geometry by a gas phase electron diffraction 



study [11]. The only fluorocarbon investigated in the supercritical regime is tetranuoromethane[12], 
whereas deuterated trifluoromethane has been investigated in the liquid regime [13, 14]. In both cases 
the total neutron structure factor has been determined. References to the simulation studies performed 
for HCF 3 can be found in the recent work by Hloucha et al [15] beginning with the early work of B6hm 
et al [16]. 

The method used here to deduce molecular orientations is based on the potential inversion scheme 
of Levesque, Weis and Reatto [17]. From the result of hard X-ray and neutron diffraction experiments 
with isotopic substitution (NDIS) a site-site potential is deduced, which in turn can be used in a NVT- 
Monte-Carlo simulation to obtain the orientational correlations. The determination of the potential is 
facilitated by the ease with which the density of these systems can be varied, both having a critical point 
at about room temperature. 

2. Experimental 

We have investigated the structure of the fluid fluorocarbons in a range of pressures (28-100 bar) and 
temperatures (298-333 K) around the critical points of 58 bar, 317.8 K for H 3 CF and 48.3 bar, 299.1 K for 
HCF3 [18]. It was aimed to combine the information of a NDIS experiment [19, 20] and a hard X-ray 
diffraction experiment [21]. The neutron and hard X-ray experiments were both carried out using the 
same sample environment (Fig. 1) built especially for this experiment. The mechanical requirements 
of the pressure cell are moderate, and aluminum is a very suitable material for the sample container. 
Aluminum has a quite low scattering power for both neutrons and X-rays, only few powder lines due 
to its cubic structure and shows only little activation in a neutron beam. The samples, DCF3 (98% D, 
Cambridge isotopes), HCF 3 , a 1:1 mixture HCF3/DCF3 and H 3 CF (all Linde technical gases) can be 
condensed into the sample container through immersion in liquid nitrogen. The cell is then mounted inside 
a vacuum tank of the neutron or hard X-ray diffractometer. The sample pressure inside the mounted cell 
can be varied via an inert gas line, separated by a steel bellow from the sample. The temperature can 
be varied with a small heater at the bottom of the cell. The price of the deuterated gases requires the 
reduction of dead volumes inside the cell: The cell can be filled with ~ 4 g sample. 

Figure 1. Schematic drawing and picture of the pressure set-up 

A: sample container, B: steel bellow, C: pressure sensor, D: heater, E: temperature sensor, F: vent for 
the pressurizing medium, G: sample vent. 

The neutron diffraction measurements were performed at the diffractometer D4b[22] at the ILL 
reactor source in Grenoble, using a wavelength of 0.7501 A. Representative examples of the raw data 
are shown in Fig. 2, showing that the aluminum cell is clearly a viable alternative to the more usual 
vanadium and titanium-zirconium cells. The hard X-ray diffraction was performed at the high energy 
bcamline BW5 at the DORIS storage ring at HASYLAB, Hamburg in its set-up for liquid and amorphous 
substances [23], using a wavelength of 0.1282 A 

3. Data analysis 

The data were corrected for systematic effects like detector dead time, absorption, container scattering, 
multiple and incoherent scattering, using the procedure described in some detail in [24], and then 
normalized. The differential cross sections are expressed in terms of the scattering functions S^(Q) 
and (Q) for the hard X-ray and the neutron cases 

S (n ){Q) = {dn) N ^ +1 (1) 

S {X) (Q) = i(Q) + 1 = [dn) 1 ; l ^ lh + 1 (2) 




Figure 2. Raw data of the neutron diffraction experiment on trifluoromethane at D4b. . 

A-C: Sample + cell scattering intensity of HCF3, MCF3 (the H/D mixture) and DCF3 respectively at 

liquid-like densities (30 °C, 100 bar), D: scattering intensity of C after subtraction of the cell scattering, 

E: Sample + cell scattering intensity of DCF3 at gas-like densities (60 °C, 32 bar), F: empty aluminum 

container. 



where (4^) is the coherent differential cross section, bi the coherent scattering lengths [25], fa the X- 
ray form factors in the independent atom approximation [26], a e i the scattering cross section of the free 
electron, Vi the stoichiometric coefficient of the atom i, and where the sums are extending over the number 
of distinct atoms in the molecule N uc , the subscript uc referring to the unit of composition, the molecule. 
In Fig. 3 the density dependence of the X-ray structure function of H3CF is shown. Beyond Q ~ 2.5 A -1 
the interference scattering intensity is dominated by the intramolecular contributions. Fitting of the 
Debye equation in the range [4 < Q < Q ma x]- 

with rij teq the equilibrium distance of the atoms i and j within the molecule and 7^ the displacement 
parameter, leads to tqf = 1.416(8) A and 7cf = 0.050(18) A independent of density. Likewise, the 
molecular parameters of trifluoromethane were determined and are shown in Table 1. For trifluoromethane 
the molecular parameters are also independent of the density and in excellent agreement with the gas 
phase values from ref. [11]. For the remainder of the article only the intermolecular contributions to the 
structure are considered. 

The intermolecular scattering contribution is related to the weighted intermolecular pair distribution 
functions by a Fourier-sine transformation: 

Q ■ (S (n) - 1) sin(Qr)rfQ (4) 
Q-i(Q)sm(Qr)dQ . (5) 
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Figure 3. X-ray structure function i(Q) of H3CF at various densities. 

A: puc = 17.6 ■ 1CT 3 A" 3 , B: p uc = 13.0 ■ 1CT 3 A" 3 , C: p U c = 10.4 • 10~ 3 A~ 3 , D: p U c = 7.3 ■ 10~ 3 A~ 3 , 
E: p uc = 2.6 ■ 10~ 3 A" 3 , F: Fit of equation 3 to the data with r CP = 1.416(8) A and 7 CF = 0.050(18) A. 

where p uc is the density per unit of composition (molecule), and and are weighted sums of the 
partial (site-site) pair distribution functions (PPDF): 

9 W =2^ w H9ij with Wij = J ' (6) 

ij y2 ^ % 1 l> 

and 

g M=^FT[ Wij (Q)]® 9ij with Wij (Q)= ^®^ Q) , (7) 

where FT() is the Fourier sine transformation and ® the convolution operation. 

Trifluoromethane has six PPDF, and four independent diffraction experiments were carried out. 
Consequently, assuming the independence of the structure from the isotopic composition, a set of one 
PPDF (HH) and three independent composite partial pair distribution functions (CPPDF) can be 
isolated, each of the CPPDF is the weighted sum of two PPDF. The three CPPDF are dominated 
by the fluorine PPDF, FF, FC and FH, while the carbon PPDF, CC and HC do contribute only very 
little. Alternatively the total pair distribution function can be split into HH, HX and XX CPPDF, with 
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Table 1. Molecular structure of HCF3. 

All distances and displacement parameter in A, all angles in degree. 

1 The angles are not refined directly, but determined from the maxima of the distance distributions. 

2 Mean value of intramolecular distances and bond angles non equivalent in the crystal. 



X either C or F. This is the same separation as that used in the case of the hydrogcnhalides, with X=C1 
in that case, thus enabling a proper comparison of our results with measurements of the structure of fluid 
HC1. All pair distribution functions are defined such that lim^oo g(r) = 1. 



4. Potential inversion 



In order to generate a three dimensional picture of the structure from the pair distribution functions, the 
potential inversion scheme of Levesque, Weiss and Reatto [17] (LWR-scheme) was applied. The idea of 
this method is based on the equation: 

—v(r) 



g(r) = exp 



+ g(r) - 1 - c(r) + B{r, v) 



(8) 



kT 

relating the pair distribution function and the pair potential, where v(r) is the pair potential, c(r) the 
direct correlation function and B(r,v) the bridge function. Starting with a first guess of the potential 
e.g. by neglecting the bridge function, a Monte Carlo simulation gives g(r)^ and c(r) ( - 1 -' belonging 
to and thus B(r, v^). Substituting B(r, z/" -1 )) for B(r,v) in equation 8 gives the LWR iteration 
formula 

y (n) j kT = v (n-l) + l n ( g («-l) /g( e *P)) 

+ c^-V - c (exp) - g( n ~V + g (exp) (9) 

Schommers [27] proposed a similar iteration scheme where only the logarithmic term of Eq. 9 is considered. 
Reatto et al have shown [28], that their algorithm converges - under certain conditions much - faster 
compared to the Schommers scheme. Soper's EPMC algorithm [29] uses the same iteration scheme 
as Schommers, but contrary to it is formulated for multi-element fluids. Likewise, it has been shown 
by Kahl and Kristufek [30] that the LWR scheme is applicable to polyatomic systems. The systems 
investigated here are even a step more complicated than the systems Kahl and Kristufek used, as the 
sites are connected by covalent bonds and the PPDF are not complete. Thus the HC and the CC site-site 
potentials were kept constant as hard sphere potentials. The method has already been tested under these 
conditions and compared to the results of Reverse Monte-Carlo (RMC) [31] simulations. A short account 
of the comparison between the potential inversion and the RMC method has been given elsewhere [32]. 
In the following only the results of the potential inversion method are given. 

One test of the potential inversion scheme is the correct or incorrect prediction of the pair distribution 
function at several different state points from a pair potential determined at a specific state point. This is 
demonstrated in Fig. 4. This procedure also tests whether the intermolecular potential can be described 
as an effective two-body potential. 



5. Results and discussion 

In Fig. 4a the density dependence of g^xx °f trifluoromethane is shown. The density dependence of the 
other CPPDF is similar. At the higher densities <?xir snows a typical liquid-like behavior and several 
maxima and minima. The maxima at larger distances die out when lowering the density, while the height 
of the first maximum increases. This is opposite to the behavior of fluid hydrogen chloride [34] or water 
[35] where the height of the main maximum decreases with density. The position of the main maximum 
remains almost unchanged. 

The HCF 3 molecular potential has been determined as a site-site potential, whereas for H 3 CF 
only the X-ray weighted pair distribution function g(r)^ was determined, and thus an independent 
determination of H 3 CF site-site potentials was not possible. The question arises as to whether the HCF 3 
site-site potentials can also be used to describe the H 3 CF structure, i.e. if these potentials have a general 
applicability to all fluorocarbons. Figure 5 indicates that the site-site potentials are indeed transferable 
to a different molecular species. 




Figure 4. a)Example of the density dependence of the intermolecular part of the g^x) CPPDF of 

HCF3. The experimental g^x ls compared with the prediction of a potential model derived at some 
specific state point (f:p uc =1.15 10 — 3 A~ 3 ). 

b) Direct comparison of the experimental g^xx a * decreasing density, (same state points as in a) 
The densities in the figure correspond to the following experimental conditions (from high to low density) 
: a: p=100 bar, T=298K, p uc = 8.8 1CT 3 A~ 3 , b: p=100 bar, T=333K, p uc = 5.8 1CT 3 A~ 3 , c: p=80 
bar, T=333K, p uc = 3.8 10~ 3 A™ 3 , d: p=60 bar, T=333K, p uc = 2.5 10~ 3 A" 3 , e: p=50 bar, T=333K, 
p„ c =1.6 1(T 3 A" 3 , f: p=40 bar, T=333K, p uc =l.l5 1CT 3 A" 3 , g: p=32 bar, T=333K, Pllc =0.84 1(T 3 

A- 3 . 
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Figure 5. Comparison of g^(r) for H3CF from the experiment with the result from a Monte-Carlo 
simulation using the site-site potentials derived for HCF3. 

The symbols correspond to the experimental g^(r), the separation of the symbols in x-direction 
corresponds to the experimental resolution, the solid line to the simulation. The g( x )(r) correspond 
to different densities shown in the figure (state point B, D and E in Fig. 3). 

The aim of the present work was to determine the influence of the molecular dipole on the 
orientational correlations between the molecules. Fig. 6 compares the HH, HX and XX pair distribution 
functions of fluid HC1 [34], HCF 3 and H 3 CF. The three PPDF of HC1 are quite structured and dissimilar 
while the corresponding CPPDF of both HCF 3 and H 3 CF are much less structured and very similar. 
This behavior is an indication that there will be no strong preference for particular orientations in the 
fluorocarbons. The most remarkable difference can be seen in the XX-(C)PPDF: gcici m HC1 has the 
highest maximum. 

Fig. 7 quantifies this qualitative statement and shows P(rcoM, cos(0)), the relative probability of 
finding a second molecule at the center of mass distance tcom in an orientation cos(0), where 9 is the 
angle between the molecular dipoles. This figure is to be compared to figures 5 and 7 of Ref. [34] and 
to figure 9 of Ref. [12]. While the first work determines the orientational correlation in fluid HC1 via 
the EPMC formalism (truncated version of Eq. 9) and finds pronounced orientational correlations, the 
second is a Reverse Monte-Carlo study using a total neutron structure factor measurement of fluid CF4 
as input and finds a P(rcoM,9) very similar to fig. 7, structured only at very short distances. 

Averaging the mean < cos(#) > r via 

/>oc 

g K = 1 + / N(r) < cos(0) > r dr (10) 
Jo 

with N(r) = 4:irp uc r 2 gcoM, to yield a Kirkwood g-factor leads to values very close to one (0.995 on 
average) in the entire density range investigated for both HCF 3 and H 3 CF. This is in agreement with 
advanced theories of the dielectric properties of these materials [33] . 

The largely simplified model illustrated in Fig. 8a can help to understand this behavior. At a distance 
of 5 A, the maximum of the pair distribution functions of HCF 3 , the energy difference between parallel 
and antiparallel alignment of point dipoles of 1.65 -3.336 • 10~ 30 Cm located at the center of mass is about 
2kT for the polar positions and about lkT for the equatorial positions. In a two level system this would 
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Figure 6. Comparison of the HH, HX and XX CPPDF of HCF 3 and H 3 CF with the PPDF of HC1[34]. 
The CPPDF for H 3 CF are simulation results 





Figure 7. a: Cosine distribution of the angle 9 between the molecular dipoles of HCF3 at p=8.S 
10-3A- 3 , 

b: Mean cosine of the angle versus the center of mass distance at three different densities. 



lead to seven times more parallel than antiparallcl orientation in polar positions and three times more 
antiparallel orientation in equatorial positions. But the equatorial region is four times larger, leading to 
an almost complete cancelation of parallel and antiparallel orientations. 

For HC1 a preference for a polar arrangement of molecules has been found [34] (molecules directly 
'below' or 'on top' of each other referred to the direction of the dipole). For these molecules in polar 
arrangement a strong preference for parallel orientations has been found - up to 17 times more parallel 
orientated dipoles than expected from a random distribution. The cosine distribution to be expected for 
point dipoles in a polar arrangement is: 

exp[fc(r) cos(0)] 



P[r, cos(0)] 



exp[fc(r) cos(9)]d cos 9 



(11) 




Figure 8. a) Simplified model of the relative orientation of dipolar molecules. Explanation in the text, 
b) Cosine distribution and theta distribution to be expected from two point dipoles in polar arrangement 
for p = 1.65 ■ 3.336 ■ 1CT 30 Cm, r ma x=5A (HCF 3 , solid line) and p = 1.07 ■ 3.336 ■ 1CT 30 Cm, r m ax = 
3.6 A (HC1, broken line) 

with k(r) = 2 • p 2 /Airer 3, . This function is shown for r — r max , with r max the first maximum of g(rcoM) 
in Fig. 8b. With p = 1.07 • 3.336 • 10" 30 Cm [36] and r max = 3.6 A thus only 3.5 time more parallel 
oriented dipoles as in a random distribution should be found. From a simple point dipole model a less 
pronounced preference of parallel orientation as found by [34] is predicted for HC1, while on the other 
hand from this model a stronger preference than actually found would be expected for both, HCF3 and 
H3CF, probably due to the detailed molecular geometry and specific site-site interactions. These two 
effects, the cancelation of parallel and anti-parallel orientations and specific site-site interactions lead 
finally to the average behavior shown in Fig. 7 very similar to fluid CF4. 

Hloucha et al. [15] recently published a constant NPT Monte-Carlo simulation of liquid HCF3 at 
subcritical temperatures. Their model is a rigid five site model, with a Lennard-Jones contribution, partial 
charges at the sites and point dipoles at the center of mass, with a constant and an induced contribution. 
With this model and in the dense liquid more pronounced orientational correlations between neighboring 
molecules as in this work are found. The positive peak in < cos(#) >,. at r ~ 4.5 A in Fig. 7 is higher by 
a factor of two and the negative peak at close contact is missing. Hloucha et al. observed a decreasing 
trend in the orientational order with decreasing density which will tend to level out the differences in 
orientational order in the range of densities investigated here. Although the detailed comparison is 
complicated by the difference in the range of thermodynamic parameters investigated, this is supporting 
the point of view that the orientational ordering in HCF3 is even less pronounced than predicted by a 
dipolar picture. 

The spatial arrangement of neighboring molecules is illustrated in Fig. 9 which compares the crystal 
(a) - the positional parameters are taken from [10] - and the fluid at liquid like dcnsities(b). In the crystal 
each trifluoromethane molecule is surrounded by twelve neighboring molecules at nearly the same distance. 
Among these, two are oriented parallel and two antiparallel, the remaining two times four molecules in 
two different T-orientations, that are orientations where the dipole moments are perpendicular to each 
other. Fig. 9b is a snap-shot of a simulation at p uc = 0.0088A~ 3 . Again, the twelve next-neighbors are 
shown. The orientation of these molecules has been grouped into four classes having a cos(#) between -1.0 
and -0.5, -0.5 and 0.0, 0.0 and 0.5 and 0.5 and 1.0, respectively, with the angle between the molecular 
dipoles. Again, quasi T-orientations occur more often. In the fluid this snap-shot is only representative 



of course, the ensemble of structures leads to Fig. 7. 



Figure 9. Relative orientation of next-neighbor molecules in crystalline HCF3 (a) and dense fluid HCF3 
(b). 

The relative orientation to the central molecule is indicated with a color code. In the crystal: Red (R): 
parallel, green (G) antiparallel, blue (B) and black (b): different T-configurations. In the fluid: red 
1.0 > cos(6») > 0.5, blue 0.5 > cos(6») > 0, green > cos(e) > -0.5, black -0.5 > cos(0) > -1.0. 



6. Conclusion 

The density dependence of orientational correlations in fluid trifluoromethane and methylfluoride has 
been studied by NVT-Monte-Carlo simulations using effective site-site two body potentials derived 
via the Levesque-Weis-Reatto inversion scheme from NDIS and hard X-ray diffraction data of fluid 
trifluoromethane. Advanced theories of the dielectric properties of these materials predict a Kirkwood-g- 
factor close to one in the entire density range investigated here. The orientational correlations found in 
the simulation are in full agreement with this prediction. Detailed comparison with the orientational 
correlations found in fluid hydrogcnchloridc and tctrafluoromethane shows that the orientational 
correlation in HCF 3 and H 3 CF are closer to tctrafluoromethane than to hydrogenchloride, although 
the interaction energy of the molecular dipoles is comparable to HC1 in these systems. This and the 
comparison with the simulation results for a dipolar model of HCF 3 suggests that in the fluorocarbons 
site specific interaction results in weaker orientational correlations than predicted by a dipolar model 
while in fluid HC1 the orientational ordering is enhanced. 
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